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Abstract 

In this paper, we describe some recent results obtained in the context of vec- 
tor subdivision schemes which possess the so-called full rank property. Such 
kind of schemes, in particular those which have an interpolatory nature, are 
connected to matrix refinable functions generating orthogonal multiresolution 
analyses for the space of vector-valued signals. Corresponding multichannel 
(matrix) wavelets can be defined and their construction in terms of a very ef- 
ficient scheme is given. Some examples illustrate the nature of these matrix 
scaling functions/wavelets. 
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1. Introduction 

This paper deals with the construction of the proper "wavelet" tools for 
the analysis of functions which arc vector-valued. Such type of functions arises 
naturally in many applications where the data to be processed are samples of 
vector-valued functions, even in several variables (e.g. electroencephalographic 
measurements, seismic waves, color images). 
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Up to now, most signal processing methods handle vector valued data com- 
ponent by component and thus ignore the possible relationship between several 
of the components. In reality, on the other hand, such measurements are usually 
related and thus should be processed as "complete" vectors. While it is obvious 
that for vector valued signals a standard scalar wavelet analysis does not take 
into account the correlations among components, even an approach based on 
the well known multiwavelets, introduced in [7] , |12j and widely studied (see [5] 
and the exhaustive references therein), is only apparently justified by the "vec- 
torization" step required in order to apply them. Indeed, multiwavelets still 
generate multiresolution analyses for the space of scalar-valued functions. 

In [T5J, [T3] matrix wavelets have been proposed for the analysis of ma- 
trix valued signals while in pQ the concept of multichannel wavelet (MCW) 
and multichannel multiresolution analysis (MCMRA) have been introduced for 
processing vector-valued data. Actually, the two concepts look very similar: 
the crucial point of both is the existence of a matrix refinable function, which 
satisfies the so-called full rank condition. However, on one hand, the matrix 
analysis scheme proposed in [T3J, [H] coincides with the applications of a fixed 
number of MCW schemes, on the other hand, only a very special characteriza- 
tion of filters is proposed, which gives rise to interpolatory orthonormal matrix 
refinable functions. Thanks to this special structure, the corresponding wavelet 
filters are easily found as in the scalar situation. Nevertheless, in [T3], [T3] the 
authors do not provide either any constructive strategy nor any example. Some 
examples of 2 x 2 matrix wavelets are derived in |15) , but they strongly rely on 
the sufficient condition for convergence described in [15] . and belong again to 
a very special class. In passing, we also mention reference [6] since, in spite of 
the title, it actually contains examples of wavelet filters specifically designed for 
vector-fields. 

In some recent papers [2J, [3J, [1] we investigated full rank interpolatory 
schemes and showed their connection to matrix refinable functions and multi- 
channel wavelets. In particular we proved that there is a connection between 
cardinal and orthogonal matrix refinable functions which is based on the spectral 
factorization of the (positive definite) symbol related to the cardinal function, 
thus, as a by-product, giving a concrete way to obtain orthogonal matrix scaling 
functions. The corresponding multichannel wavelets, whose existence is proved 
in pQ , can be found in terms of a symbol which satisfies matrix quadrature mirror 
filter equations. 

In this paper, we show how the solution of matrix quadrature mirror fil- 
ter equations can be found. In particular we give an efficient and construc- 
tive scheme which makes use again of spectral factorization techniques and of 
a matrix completion algorithm based on the resolution of generalized Bezout 
identities. 

The paper is organized as follows. In Section 2 all the needed notation and 
definitions are set. Some introductory results on full rank vector subdivision 
schemes are also presented and their connection to orthogonal full rank refinable 
functions is shown. In Section 3, the concepts of MCMRA and MCW are revised 
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while their orthogonal counterparts are considered next in Section 4. Section 5 
deals with the construction of multichannel wavelets and an explicit algorithm 
for such construction is presented. In Section 6, starting from some interpolatory 
full rank positive definite symbols, the associated multichannel wavelets with 
two and three channels are constructed and shown. Some conclusions are derived 
in Section 7. 



2. Notation and basic facts on vector subdivision schemes 

For r, s e N we write a matrix A <G W xs as A — (ajk ■ j = 1, • • • , r, k = 
1, . . . , s) and denote by £ rxs (Z) the space of all r x s-matrix valued bi-infinite 
sequences, A = (Aj e W xs : j e Z). For notational simplicity we write £ r (Z) 
for £ rxl (Z) and £(Z) for i 1 (Z) and denote vector sequences by lowercase letters. 

By ^2 XS (^) we denote the space of all sequences which have finite 2 norm 
defined as 




(1) 

where |-| 2 denotes the 2 operator norm for r x s matrices. 

Moreover, L r 2 xs (M) will denote the Banach space of all r x s-matrix valued 
functions on R with components in and norm 

\ 1/2 

EE/i^wi 2 ^ ■ (2) 
K k=i 3=1 j 

For a matrix function and a matrix sequence of suitable sizes we introduce 
the convolution "*" defined as, 

(F*B) :=Y J F{--k)B k . 

fcez 

For two matrix functions F, G € L r 2 xs (R), their inner product is defined as 



(G,F) = [ G H (x)F(x)dx. 



and we have 

||F|| 2 = (trace (F,F))* . 

Let now A be a matrix sequence A = (Aj e R rxr : j e Z) E f rxr (Z). Let c 
be any bi-infinite vector sequence c = (cj e R r : j G Z). The vector subdivision 
operator S^, based on the matrix mask A, acts on c by means of 

Sac = I E A 3-2k c fc : j e Z J . 
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A vector subdivision scheme consists of iterative applications of a vector 
subdivision operator starting from an initial vector sequence c € £ r (Z), namely: 

c° :=c 

c n := S^c"" 1 = S^c° n > 1. 

It is L2- convergent if there exists a matrix-valued junction F £ L^* 1 *^) such 
that 

lim 2-™/ 2 ||/i"(F) - S^I|| 2 = 0, (3) 

n— »oo 

where the mean value operator at level n £ N is defined as 

(fj, n (F))j := T [ F(t)dt, j e Z. 

J2-"(j + [0,l]) 

The matrix-valued function F associated in such a way with a convergent sub- 
division scheme is called the basic limit function and it is refinable with respect 
to A, that is 

The symbol of a subdivision scheme Sa or of the mask A is defined as: 
A(z)=J2A jZ j , zeC\{0}, 

and the subsymbols, A £ (z) = X^ez A s+2j z° > z G C\{0}, £ € {0, 1}, are related 
to the symbol by 

A(z) - Ao(z 2 ) + zA!(z 2 ), zeC\{0}. (4) 
Next, we define the two r x r matrices 

A e := A e (l)=^A e+2i , £€{0,1} 

jez 

and their joint 1-eigenspace 

S A ■= {ye« r : A a y = A x y = y} . 
The rank of the mask A or of the subdivision scheme Sa is the number 

R(A) := dim£ A 

satisfying 1 < R(A) < r for convergent schemes. In particular, is said to be 
of full rank if R(A) = r. 

The following proposition gives equivalent characterizations of full rank schemes 

m 

Proposition 1. Let Sa be a vector subdivision scheme with symbol A(z). Then 
the following statements are equivalent: 
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i) Sa is of full rank; 

ii) the symbol satisfies .4.(1) = 21 and A(—l) = 0, or, equivalently, there 
exists a matrix mask B £ £ rxr (Z) such that A(z) = (z + l)B(z), with 
B(l) = I (scalar-like factorization); 

iii) the scheme preserves all vector constant data, i.e whenever Cj = d, j € Z, 
for some d £ W , then also (S^c)j = d, j £ Z; 

iv) the basic limit function F is a partition of the identity, i.e. F(- — j) = I. 

A full rank vector subdivision scheme is interpolatory if it is characterized 
by the property 

{S A c) 2j = Cj, j £ Z, 

or, equivalently, 

A 2j = S j0 I, j £ Z. 

Its associated basic limit function, whenever it exists, turns out to be cardinal, 
i.e., 

F(j)=6 j0 I, j£l. 

The connection between interpolatory and full rank subdivision schemes is ex- 
pressed by the following result [2]: 

Proposition 2. An interpolatory subdivision scheme with symbol A(z) is of 
full rank if and only if A(l) = 21. 

3. Multichannel multiresolution analysis and multichannel wavelets 

As in the scalar situation, a multichannel multiresolution analysis (MCMRA) 
in the space L£ (M) of square integrable vector valued functions can be defined 
b a nested sequence ■ • • C VL_i C Vo C V% C • • • of closed subspaces of L\ (R) 
with the properties that 

1. they are shift invariant 

h£Vj h(- + k) G Vj, k£Z; 

2. they are scaled versions of each other 

h£V <S=^> h(2 j -)£Vf, 

3. f)V 3 ={[0,0 ...Of}, \Jv J =L r 3 (R); 

4. the space Vq is generated by the integer translates of r function vectors, 
that is there exist /' = (/|, . . . , /*) £ L\ (R), i = 1, . . . ,r, such that 

i) Vq = span {/*(• — k) : k € Z, i = 1, . . . , r}, 
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ii) there exist two constants K\ < K 2 such that 

/ r \ V2 



MEE 

\j=i feez 



ft 



< 



EE4/ J '(--fc) 



1=1 



1/2 



<MEE 



(5) 

for any r scalar sequences c 1 — {c k } G £2(^)1 i = 1, • • • , r. 

Now, let h = (hi, . . . , h r ) T G Vo C L^M.). Then there exist r scalar se- 
quences c 1 = (c\ G R, : fc G Z) G ^(Z), i = 1, . . . ,r, such that 

^ = E( c ^ 1 (-- fc ) + ---+ c ^ r (-- fc ))' 



feez 



that is 







/ 








f /r(- 


-fc) ^ 


\ 




= E 


4 




- k) ) + 


• • • + 4 








V K J 


fcez 




{ fr(- 


-k) ) 




V /r r (" 


-fc) y 


/ 



/ //(--fc) ••• /f(--fc) \ / 4 \ 



= E 

feGZ 



V fr(-k) ■ ■ ■ /;(■ -fc) 

F(-fc) 



4 / 



Thus, any ft, <G Vq can be written as h = F * c — E -^(' — fc) c fc, where F G 



L r 2 xr { 
have 



and c = (c fe G R r : fc G Z) G ^(Z). Since f 1 , . . . , f G V C 14, we 

r = E^(2--fc)4, i = l,...,r, 

fcez 

for some vector sequences a 4 = (a\ G R r , : fc G Z) G ^2(^)1 « = 1, ■ • ■ , Thus, 
the matrix-valued function 

F = (/ i |...ir) G^ xr (R) 

satisfies the matrix refinement equation 

F = Y J F(2--k)A k 



(6) 



where *4 = = (a\ \ . . . \a r k ) G R rxr : fc G Z) is the refinement mask. 

Now, since the subdivision scheme £4 applied to the initial sequence SI 
converges to a matrix function with r stable columns it follows that £4 = R r 
so that 

A £ = E^-2j =/, £G {0,1}, 



G 



i.e. A(l) — 21, A(— 1) = 0, which is the full rank property of the mask A. 

The "joint stability" condition ^ on the function vectors f\ j = 1, . . . , r, 
implies that the function F is stable in the sense that there exist two constants 
K\ < K 2 such that 

K x ||c|| 2 <||F*c|| 2 <^ 2 ||c|| 2 , ce*£(Z). 

Summarizing the described properties, the following important result holds 
true: 

Theorem 1. The subspaces Vj C L 7 2 (R), j e Z, /orm an MC'MRA if there 
exists a full rank stable matrix refinable function F € L2 Xr (ffi) such that 

V 3 = {F*c(2 J '-) : c€*£(Z)}, j€Z. 

In analogy to the scalar case, we call F matrix scaling function. 

4. Orthogonal multichannel mult iresolut ion analysis 

For many application purposes, orthogonal MCMRAs are the most interest- 
ing. They are are equivalently characterized by the following properties: 

1. the space Vo is generated by orthonormal integer translates of the function 
vectors f J ,j = l,...,r; 

2. the matrix scaling function F is orthonormal that is 

(F,F(- + k)} := [ F H {x)F{x + k)dx = S k0 I, k e Z. 

The key for the construction of orthonormal matrix scaling functions (and 
the associated multichannel wavelets) is their connection with interpolatory vec- 
tor subdivision schemes, as established in the following Theorem where the 
canonical spectral factor A{z) of a positive definite interpolatory full rank sym- 
bol C{z) is the symbol satisfying C(z) = \ A H (z)A(z) and A(l) = 21, A(-l) 
0. 

Theorem 2. [2] Let C(z) be the symbol of an interpolatory full rank vector 
subdivision scheme. Let C{z) be symmetric strictly positive definite for all z ^ 
— 1 such that \z\ = 1. Then the canonical spectral factor A{z) ofC(z) defines 
a subdivision scheme which converges in L2 Xr (M.) to an orthonormal matrix 
scaling function F. 

As to the construction of such function, in [2] a simple procedure is given 
based on a modified spectral factorization of the parahermitian matrix 2C(z), 
which takes into account the presence of some of its zeros on the unitary circle. 

We can now associate a matrix wavelet to any orthogonal MCMRA in the 
standard way. A matrix function G € V\ is called a multichannel wavelet for 
the orthogonal MCMRA if: 
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1. Wj := V j+ i eVj = {G*c(2>-) : ce F 2 (Z)} , j € Z; 

2. G is orthonormal. 

As to the existence of multichannel wavelets, we recall the following result: 

Theorem 3. f^/ Suppose that the matrix scaling function F G L2 Xr (R) has 
orthonormal integer translates. Then there exists an orthonormal wavelet G G 
L2 Xr (IR) satisfying the two-scale relation 

G = Y J F{2--j)B J (7) 

for a suitable mask B — (Bj : j £Z). 

Observe that, as in the scalar case, due to the full rank properties of F, the 
symbol of the multichannel wavelet G must possess at least one factor [z — 1). 
This is equivalent to say that the multichannel wavelet has at least one vanishing 
moment. 

It is easy to show that the symbols of an orthonormal matrix function F and 
of the corresponding multichannel wavelet G satisfy the orthogonality ( quadra- 
ture mirror filter ( QMF)) conditions 

A\z)A{z) + AH-z)A(-z) = 41, \z\ = l, 
A t (z)B(z) + A i {-z)B(-z) = 0, \z\ = l, 
B\z)B{z) + B\-z)B{-z)=AI, \z\=l, 

where A^[z) := A T (z~ 1 ) which means A^(z) :— A H (z) whenever \z\ = 1. 

The QMF can be written in a concise form by the condition U^(z)U{z) = I 
on the block matrix 

TTf . If A{z) B{z) \ , rjU . 1 / A\z) A\-z) 
U( - Z) := 2 [ A(-z) B(-z) ) ' Where U "M 2 [ B\z) B\J) 

(which is the condition of being an unitary matrix for \z\ = 1). 

Note that, given an orthogonal symbol A(z), the alternating flip trick, used 
to construct the wavelet symbol in the scalar case, does not work in this context. 
Indeed, the symbol B(z) = zA^(-z) does not verify the QMF equations unless 
A{z) and A(—z) commute. 

On the other hand, the orthogonality conditions can be written in terms of 
the subsymbols Aq(z), Ai(z) and Bq(z), B\{z). In fact, from Q we see that 
the matrix V[z) 

Viz) - ( Ao ^ )- ( 1 M U(z) 

is such that V^(z)V(z) — 21. Thus, the QMF equations take the equivalent 
form 

^(z 2 )A (z 2 ) + A\(z 2 )Ai(z 2 ) = 21, \z\ = 1, 
Al(z 2 )B (z 2 ) + A{(z 2 )B 1 (z 2 )=O, \z\ = l, (8) 
B\ (z 2 )B (z 2 ) + Bliz^B^z 2 ) = 21, \z\ = 1. 
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Suppose now that we are given an orthogonal MCMRA generated by some 
orthonormal matrix scaling function F, with G as the corresponding multichan- 
nel wavelet. Let us consider any vector- valued function h E L^M). From the 
nesting properties of the spaces {Vj} and {Wj}, it follows that the approxima- 
tion Pfh of h in the space Vg, £ E Z, can be found in terms of the following 
multichannel wavelet decomposition: 

P e h = P e - L h + Qe-ih + Q t - 2 h + ... + Qe-Lh, 

where L > and Pg-^h, Qz_jh, j = 1, . . . , L, represent the orthogonal projec- 
tions of h to the spaces Vi-L,Wg-j 7 j — 1, . . . , L, respectively. Analogously to 
the scalar case, we can derive a fast algorithm which allows to compute all the 
projections by means of a recursive scheme. In fact we have that, if Pjh E Vj 
then 

kez 

where 

c[ j) = {F{2° ■ -k),h) = J F H (2 j x - k)h{x)dx. 

To compute the vector coefficient sequence c^ -1 ^ E l 2 {1) connected to the 
representation of h in the space Vj-%, we make use of the refinement equation 
and get 

(F(2^ 1 • -k),h) = (J2 F{2{V- 1 ■ -k) - n)A n , h) 
J2^l-2 k {F(2 j ■ -n),h) = XX_2*4> k e Z- 

nGZ ngZ 

In the same way, using the wavelet vector coefficients sequence d^~ 1 ^ E 
^(Z) connected to the representation of h in the space Wj-i is obtained as 

d^- x) = (G(y- 1 --k),h)=J2B? l _ 2k c> n , kEl. 

In summary, assuming I = 0, that is Vo as the initial space of our representation, 
the vector decomposition formula up to the level L > reads as 

c^- 1] =J2 A n-2k<i, 4 j - 1} = J2B^ 2k c> n , fc€Z,j = L. 

Conversely, given the projections Pjh and Qjh, the vector coefficient sequence 
c 0+!) connected to the representation of h in the space Vj+i = Vj © Wj is 
obtained by considering that 

p J+1 h = Y. F ^ +1 -- k ) c i 3+1) (9) 

kez 



,07-1) _ 
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and, on the other hand, 

P j+1 h = Pjh + Qjh = F(V ■ ~n)c^ + ^ G{V ■ -n)t$\ 

By invoking again the refinement equation on the right-hand side of the previous 
expression, we have that 

P j+1 h = J2J2 F ( 23+1 - k)A k _ 2n c^ G{23+1 - n)B k _ 2n d^. (10) 

A comparison between the two expressions {9} and ([Tof gives the vector recon- 
struction formula: 

c fc +1) = J2 ( A k-2nc^ + B fe _ 2 „d^) , k e Z, j = —L, . . . , 0. 
For similar matrix wavelet decomposition and recontruction schemes see |14) . 

5. Multichannel wavelet construction 

The equation A\ ) {z)Aq{z) + A\(z)Ai{z) = 21 is the starting point of the 
procedure that we propose to construct the symbol B(z) of the multichannel 
wavelet. More in detail, the procedure derives the matrix Laurent polynomials 
B (z) and B t (z), and thus B(z) = B {z 2 ) + zB^z 2 ). 

In the simple situation where A (z) and A\{z) are diagonal symbols, the 
multichannel wavelet subsymbols B (z) and Bi(z) are constructed by the re- 
peated application of a scalar procedure: for each couple (ao(z))a, (ai{z))a, 
i — 1, • • • , r, the two Laurent polynomials (bo(z))u, (bi(z))u solution of the Bc- 
zout identity are derived (see Lemma [T]). 

In the general situation, the strategy is based on the existence of two matrix 
symbols Dq(z), -Di(z) such that the positive definite Hermitian matrix 

G{Z) -{MZ) D{(z)) 

satisfies 

dctG(,) = l, G»(,)G W = ( 2 J I)o(z)jD » ( , ) ( ; jDi(z) ^ i(z) ), (12) 
as proved in the following proposition. 

Propositions. Let Aq(z) and A\(z) be such that Aq(z)Aq(z) + A\ (z)Ai(z) — 
21. Suppose that there exist two matrix symbols Dq(z), Di(z) such that the 
positive definite Hermitian matrix G(z) in (11) has the properties |7 

K(z) be the spectral factor of D(z) := 2 (d q (z)D\(z) + D x (z)D\(z) 
the matrix symbols 

B {z) = 2Dl(z)(K\z)y\ B x {z) = 2D\{z)(K*[z))~\ (13) 

satisfy the two last equations in Q). 
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Proof. Since D(z) is positive definite with determinant equal to 1, its spectral 
factor K(z) exists and satisfies D(z) = K(z)K^(z). Therefore, if we substitute 
([13} into 

Al(z)Di(z) + A{(z)D[(z) = 

we get 

A\ (z 2 )B q (z 2 ) + A^B^z 2 ) = 0, \z\ = 1, 
which is the second equation in pj). Finally, since 

(K(z)y 1 D(z) (k\z))~ X =I, 



by substituting ( 13 ) into it we end up with 

B\ (z 2 )Bo(z 2 ) + B\{z 2 )B x {z 2 ) = 21, \z\ = 1 
which is the last equation in □ 

For the proof of the existence of the Hermitian matrix G(z) we refer the 
reader to |llj . Its actual construction can be carried out as explained in the 
following subsections. 

5.1. Some matrix completion results 

We start by recalling a result about Bezout identities 

Theorem 4. JiOJ / For any pair of Laurent polynomials ai(z) 02(2) the Bezout 
identity 

a x {z)bi{z) + a 2 (z)b 2 (z) = 1 (14) 

is satisfied by a pair b\{z) b 2 {z) of Laurent polynomials if and only if a\(z) and 
a 2 (z) have no common zeros. Moreover, given one particular pair of solutions 
b*(z), b 2 (z) the set of all solutions of (14-) is of the form 

bi(z) = b*(z) +p(z)a 2 (z), b 2 (z) = b* 2 (z) - p{z)a x {z) 

where p(z) is any Laurent polynomial. 

As to the common zeros of the subsymbols, the following result holds true. 

Lemma 1. The subsymbols associated with the diagonal entries of an orthogonal 
symbol A(z), that is the subsymbols (ao(z))a, {a\[z))a, i = l,---,r, have no 
common zeros. 

Proof. The result follows from the first equation in (JsJ) . In fact the existence of 
a common zero for (ao(z))a, {ai{z))u, i = 1,- • • ,r contradicts 

(a Q (z~ 1 )) li (a (z))ii + (a 1 (z~ 1 ))ii(ai(z))ii = 2, i = l,...,r, 

a relation satisfied by the diagonal symbols. □ □ 
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We continue with a "completion" result. 

Theorem 5. Let di(z), i = l,...,n, Laurent polynomials with n > 2 with no 
common zeros. Then there exists a nxn matrix Laurent polynomial P(z) whose 
first row is (ai(z), ■ ■ ■ ,a n (z)) such that 



detP(» = l, zeC\{0}. 
Proof. The proof is by induction. Let us start with n 
tor (ai(z) 1 a 2 (z)) we construct the matrix P(z) = 



with 



= 2. For the row vec 
ai(z) a 2 {z) 
-h(z) b 2 (z) 

bi(z), i — 1,2 being solutions of the Bczout identity (14). Obviously, the matrix 
P(z) is such that detP(z) = 1, z E C \ {0}. Next, assuming that the theorem 
is true for n — 1, we prove it for n. Thus, given (ai(z), • • • , a n (z)) we first con- 
struct P{z), det P(z) = 1, with the first row given by (ai(z), • • • , a n -\(z)) then 
we construct the block matrix 



P(z) 



This completes the proof. 



\ 



P(z) 



a n (z) \ 




□ 



□ 



With this result we are now able to describe a procedure for "completing" 
an m x n matrix Laurent polynomial with m > n and rank A(z) = m, by 
constructing an n x n matrix Laurent polynomial P(z), whose first m rows 
agree with those of A(z), such that 

detP(z) = l, zeC\{0}. 

The existence of such polynomial has been formerly proved in |llj . 

We start by discussing how to write a procedure providing, for a given Lau- 
rent polynomial vector a n (z) of length n > 2, a Laurent polynomial matrix 
P{z) having the first row given by a n {z) = (ai(z), • • • , a n (z)) and satisfying 
det P(z) = 1, z G C\ {0}. 

For ai(z), a 2 (z), let [bi(z), b 2 (z)]=Bezout(ai(z), a 2 (z)) be the procedure in 
Matlab-like notation providing the two Laurent polynomials b%(z), b 2 (z) solution 
of the Bezout identity. 

Next, let [P„ xn (z)]=Basic_completion(a n (z)) be the following (recursive) 
procedure that, taking the vector a(z) as input, produces the matrix P nxn (z) 
as output. 

• [Pnx7i(^)]=Basic_completion(a n (2:)) 

If n = 2 let [bi(z), 62(2)]— Bezout(ai(z), a 2 (z)) two Laurent polynomials 

ai(z) a 2 (z) 



solution of the Bezout identity. Set P2x2(^) 



M-A b 2 (z) 



else (n > 2) use [P(„_i) x („_i)(z)]=Basic_completion(a n _i(z)). 



12 



Then with e„_i := (0, • • • , 1) T e E n_1 , define 

p ._ f -P(n-l)x(n-l)(^) an(z) 

nxn\^) ■ — n 

V ^lx(n-l) e„_i 

We continue by describing the (recursive) procedure [P„ lxm ]=Completion(A„ xm ). 
This procedure, for a given matrix Laurent polynomial A nxm (z) with m > n, 
constructs the Laurent polynomial matrix P(z) whose n first rows agree with 
those of A nxm and such that det P(z) = 1, z € C \ {0}. 

• [P mX m(z)]=Completion(A„ xm (z)) 

1. Take A( n _ 1 ) xm (z) the sub-matrix of A nxm (z) made of its first (n — 1) 
rows ; 

2. If n = 2 construct [P mx _ m (z)]=Basic_completion(A lxm (z)) 
else (n > 2) construct [P mxm (z)]=Completion(j4 ( -„_ 1 ) xm (2)) ; 

3. Compute the matrix C nxm (z) := A nxm (z)P mxm (z) having block struc- 
ture 



C nxm (z) — 



I (n-l)x(n-l) 0(n-l)x(r+l) 
c lx(n-l)( z ) d lx(r+l)( z ) 



where r := m — n; 

4. Construct the (r + 1) x (r + 1) matrix 
[-D(r+i)x(r+i)(z)]=Basic_completion(d lx(r+1) (z)) ; 

5. Use -D( r+1 ) X ( r+1 )(z) to construct the mx m block matrix 





f I(n-l)x(n-l)(z) 


0( n -l)x(r+l) 


V mxm{ z ) '■ — j 


Clx(n-l)iz) 
\ Orx(n-l) 


D(r+l)x(r+l)(z) 



6. Set P mxm (z) : — V r mxm (z)J 5 mxm (z) . 

With the help of the two previous recursive procedures, we are able to sketch 
the algorithm for constructing multichannel wavelets. 

5.2. The MCW construction algorithm 

1. From A(z) e £q XT (Z) extract the sub-symbols A (z) and A 1 (z); 

2. If Ao(z) and A\(z) are diagonal symbols construct Bq(z) and B\(z) by 
the repeated application of the procedure 

[{bo(z))u, {b-L{z))u]=Bezont((ao(z))u, (ai i = 1, • • • ,r 
Then go to step 10; 

3. Otherwise construct the rx2r block matrix A (z) = ( ; -Ai(-z) ); 

4. Use [-L(-j)]=Completion(j4'(^)) to construct the matrix 

HZ) - { C (z) d(z) 



so that 

--It 

det 



A (z) Cl(z) \ 
AiW C»(*) 
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5. Compute the r x r matrix R(z) = \ \ J^(z)C\{z) + A\{z)C\{z)\ ; 

6. Compute the r x r matrices Di{z) = —B}{z)A\{z) + Ci(z), £ = 0,1; 

7. Compute the r x r positive definite matrix 

D(z) = 2 (d (z)dI(z) + D 1 {z)D{{z)) ; 

8. Compute the spectral factorization of D(z) such that 

D(z) = K(z)K i (z), 

for example using the algorithm described in [S]; 

9. Set E(z) — K~ 1 (z) and construct the r x r matrices 

B l (z) = 2D\{z)E i {z) 1 £ = 0,1; 
10. Construct the wavelet symbol 

B(z)^B (z 2 ) + zB 1 (z 2 ). 
Two remarks are worth to be made: 

• For a proof of the correctness of all steps we refer to [TTJ Section 6] . Never- 
theless, since |llj deals with rank-1 matrix functions, the above construc- 
tion has to be considered as an extention to the full rank case; 

• The algorithm produces correct results even in case of diagonal subsym- 
bols Aq(z) and A\(z). Though, to end up with exactly diagonal subsym- 
bols B (z) and Bj(z) we may need to multiply the result of the MCW 
algorithm by a suitable permutation matrix J. 

6. Numerical Examples 

The aim of this section is to consider two orthonormal full rank refinable 
matrix functions with r = 2 and r = 3, respectively and to use the MCW 
algorithm described in the subsection |5.2| to construct the corresponding mul- 
tichannel wavelets. In particular, while for the case r = 2 the positive definite 
full rank symbol that we start with is the one considered in [2], the case r = 3 
is here derived from scratch. 

6.1. A two-channel example 

Let C(z) be the symbol of the interpolatory full rank vector subdivision 
scheme first given in [2J 
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Figure 1: Orthonormal matrix scaling function with r = 2 



with 

Cn(z) = — - — — — — (linear B-spline symbol), 

t ^ 1 (^ 2 -4z + l)(z + l) 4 . 

C22{z) = ^ (4-pomt scheme symbol), 

16 z 6 

c 12 (z) = Az(z 2 -1) 3 , 

satisfying C(l) = 21, C(z) + C(—z) = 21. In order to assure positive defi- 
niteness, the parameter A is taken in the interval (-^v3, mv3)- Fixing, for 
example, A = ^j, the symbol A(z) of the associated full rank orthonormal 
rcfinablc function constructed with the algorithm given in [5] has elements 

011(2) = 1.081604742 + 0.7776021479 z + 0.2165957837z _1 - 0.08257171989z~ 2 
+0.004835070286z" 3 + 0.0009670018466z -4 + 0.0009670219776z -5 , 

012(2) = -0.001208777472 + 0.001208777472 z 2 + 0.001208777472 z - 0.001208777472z 

021(2) = -0.02117594497- 0.005136018632 z- 0.07331687185z _1 + 0.2747671662z~ 2 
-0.0679559177z~ 3 - 0.2535912212z~ 4 + 0.1464088082z~ 5 , 

a 22 (z) = 0.3169890016 + 0.6830110222 z 2 + 1.183011034a - 0.1830110103z _1 . 

The plot of the corresponding matrix scaling function is given in Fig. [T] 

Using the MCW Algorithm we find that the symbol B(z) of the correspond- 
ing orthonormal two-channel wavelet has the coefficients Bk, k = —2,..., 7, 
given in Table [I] The plot of the wavelet is represented in Fig. [2] 
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Figure 2: Orthornormal 2-channcl wavelet 



Table 1: Coefficients of the 2- 

Bk 



channel orthonormal wavelet 

k B k 



-2 



-1 



0.8140199 
-0.0014406 
-1.1322992 











0.1915220 0.0032939 

0.0037742 -0.0000058 , 

0.1360351 -0.0002724 
-0.0018743 0.6524543 

-0.0065259 -0.0133465 

-0.0003995 -1.1300526 



-0.0044401 0.0519219 
-0.0000359 0.331059 
0.0009821 0.2546764 
-0.0000192 0.1253084 
0.0007061 0.1209272 
-0.0000047 0.0153323 
/ -0.2427089 \ 
\ 0.0047518 ) 
( -0.1744916 \ 
I 0.0011525 / 



6.2. A three-channel example 

We now construct a 3 x 3 positive definite parahermitian interpolatory sym- 



bol 



cu{z) c 12 {z) c 13 (z) 
C{z) = ( ci 2 (z" 1 ) c 22 (z) c 23 {z) 

C13(Z _1 ) C23(Z~ 1 ) C 33 {z) 
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Table 2: Coefficients of the 3-channcl orthonormal matrix scaling function 



-3 






0.04 






0.04 - 

( -0.0022120 
0.0000731 
-0.12 

/ 0.0087982 
0.1481828 
0.12 

/ -0.0003396 
-0.2568731 
0.12 

/ -0.0193938 
-0.0672249 
-0.12 



V 



V 





03125 

\ 

0.03125 

-0.0017281 \ 

0.0000571 

-0.09375 

0.0060117 \ 

-0.0017879 

0.09375 

0.0004183 \ 

0.0059962 

0.09375 

-0.0240319 \ 

-0.0009911 

-0.09375 



-2 



-1 



-0.0715410 
0.2697914 
-0.04 

( 0.2874359 
0.0756559 
0.04 



1.0740925 
-0.0129914 




0.7231598 
-0.0053020 




0.0216501 
-0.0206268 
-0.03125 

0.0253824 
-0.1655344 
0.03125 

-0.0422890 
0.3501563 




-0.0073621 
1.1683135 


/ 0.0219488 \ 
0.6644173 




> 


°y 

o\ 



°y 

o\ 



o\ 



1 1 



satisfying C(l) = 21, C(z) + C(—z) = 21. It means that, as in the previous 
example, the diagonal elements should be the symbols of interpolatory scalar 
schemes, with factor (z + 1) of order mi, 1712,1713, respectively. The off-diagonals 
should certainly contain a (z 2 — 1) factor and must satisfy cy(z) = — cy(— 2) (see 
ID). We take on the diagonal the Deslaurier-Dubuc filters with m 1 — 2, m 2 = 
4, m 3 = 2, that is: 

1 (z + l) 2 1 (z 2 -4z + l)(z + l) 4 

c n( z ) = 2 — " — ' C22 ( z ) = _ Ye ' ° 33 ^ = Cll ^ z >- 

Since (z + l) 2 is a common factor, repeating the consideration done in [2], 
we require the following symbols on the off-diagonal 

c 12 (z) = Ai z (z 2 - l) 3 , c 13 (z) = X 2 z (z 2 - l) 4 , c 23 (z) - A 3 z (z 2 - l) 4 . 
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Table 3: Coefficients of the 3-charmcl orthonormal 

B k 



matrix wavelet 
k 



V 

/ 0.0297628 
1.1269644 
i 0.0367180 

/ -0.0940031 
-0.3500526 
v -0.0367180 

/ -0.2583832 
-0.1287475 
v -0.0096495 

/ -0.0908757 
-0.0045318 
0.0096495 

/ 0.2512676 
0.0029602 
l 0.0001092 



-0.0031659 

0.0001046 





-0.6416072 





V 



-0.7705220 
0.0254539 
0.0169580 

1.1447719 
-0.0169113 
-0.0169580 

-0.2560664 
-0.0115762 
-0.0033517 

-0.1606346 
-0.0170932 
0.0033517 

0.0347952 
-0.0035235 
0.0004847 



/ 







q\ 



°) 

-0.0001305 \ 


-0.9720037 

\ 

-0.0264419 
0.9720037 

-0.0751903 \ 
0.0489688 
-0.0409056 

0.0018325 \ 
-0.0567853 
0.0409056 

0.1508696 \ 
-0.0337223 
0.0160288 
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/ 

V 

/ 

V 

/ 

V 



0.1684248 
-0.0017438 
-0.0001092 

-0.0017816 
-0.0009269 


-0.0012167 
-0.0015390 


-0.0000291 
-0.0008803 


/o 










0.0212803 
0.0229928 
-0.0004847 

-0.0078392 
0.0105339 


-0.0056664 
-0.0062807 


-0.0001188 
-0.0035956 


0.0295229 
0.0373432 


0.0007057 
0.0213611 




0.0595696 
0.1406681 
-0.0160288 

-0.1305600 \ 
0.0451634 
-0.0026508 j 

-0.0773112 \ 
-0.1222042 
0.0026508 j 

0.0406917 \ 
-0.0543552 




The following parameter choice 

Ai = 1/20, A 2 = 1/50, A 3 = 1/64 

gives positive dehniteness. In this case, the canonical spectral factor A(z), 
symbol of the orthonormal refinable function F ', has the coefficients given in 
Table [2j The plot of the corresponding scaling function is represented in Fig. [3} 
The MCW algorithm applied to this symbol, gives the orthonormal wavelet 
whose symbol coefficients are given in Table [3] and whose plot is represented in 
Fig. 
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7. Conclusions 

This paper discusses a way to construct orthonormal multichannel wavelets 
from convergent full rank orthogonal symbols. As to our knowledge, the explicit 
algorithm for multichannel wavelet construction (the MCW construction algo- 
rithm) here derived, based on the strategy for rank-1 filters given in is the 
first algorithm proposed so far in the literature. 

As already pointed out, multichannel wavelets provide an effective tool for 
the analysis of multichannel signals, that is vector-valued signals whose com- 
ponents come from different sources with possible intrinsic correlations, for ex- 
ample seismic waves, brain activity (EEG/MEG) data, financial time series, 
color images. Many multichannel signals exhibit a high correlation which can 
be revealed and exploited by a MCW analysis, with filters suitably tailored to 
the specific data and application. Future investigations include the construction 
of "data- adapted" MCW bases and their application in many problems where 
vector-valued signals have to be processed for compression, denoising, etc. 
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